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Complex extensión of quantum mechanics and the discovery of pseudo-unitarily invariant random matrix 
theory has set the stage for a number of applications of these concepts in physics. We briefly review the 
basic ideas and present applications to problems in statistical mechanics where new results have become 
possible. We have found it important to mention the precise directions where advances could be made if 
further results become available. 
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1 Introduction 

The postúlate of quantum mechanics which makes a contact with measurements is the existence of real 
eigenvalues of a linear operator corresponding to a physical observable JT|. A complex extensión of quan- 
tum mechanics has been realized where certain special non-Hermitian operators are shown to possess 
real eigenvalues |2] [3] |4] [5] . Fundamental concepts like uncertainty relations in quantum mechanics |6] 
and Poincaré-Cartan integral invariant [7 1 have been generalized. A number of examples and formal results 
have been obtained [8|, matrix representations, pseudo-unitary symmetry and group structure has also been 
established [9|. On this basis, random matrix theory was presented 191 ÍTÜ1 . We are thus given a powerful 
combination of "non-hermiticity" and "random matrices" to model and describe various physical phenom- 
ena. The relevance of non-hermiticity has been well-known in many-body theory ifTTI . non-equilibrium 
statistical mechanics lTT2in~3l . open quantum systems [ 14], and number theory 1ÍT51 . 

Non-hermiticity of an operator is usually accompanied by a physical interpretation in terms of dissipa- 
tion. At a rigorous level, the breakdown of Hermiticity naturally occurs when a part of the system is not 
taken into account, as it corresponds to certain irrelevant variables. Mathematically, this was fortified as the 
Feshbach projection operator technique. This method has been very successfully employed, particularly in 
nuclear physics IÍTTI . This development has led to an understanding of many aspects of nuclear structure 
and reactions via what is known as "shell model embedded in continuum". The relevance of this is over 
the entire many body theory. To understand this, we just need to recall the first nontrivial instance where 
a localized state is coupled to a non-interacting fermionic system via one-body operators. The answer 
found there goes on to build our understanding of the Anderson model of a magnetic impurity residing in 
a metallic host iflól . 

From the first treatment of Brownian motion 11771 . random walks play a very important role in topics 
as diverse as polymer physics IÍT81 . to locomotion of bacteria IÍ191 |2ÜI . Systems out of equilibrium are 
described by a master equation with a non-hermitian operator, exemplified by the transfer matrix of the 
random walk problem. We will review the role played by the recent advances in RMT to this subject. 

RMT is also intimately connected to exactly solvable models, recognized by Sutherland lETI . The 
generalization of the projection method employed to establish this connection was accomplished and new 
model was presented Í22l for the pseudo-Hermitian case. Of great importance is the random Ising model - 
we shall argüe here that advances in RMT are needed to solve this problem. 
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2 RMT for 2x2 pseudo-Hermitian matrices 

In any physical experiment what we measure is a (real) eigenvalue of some operator belonging to a normed 
vector space. Let us consider a linear wave theory (quantum mechanics) where a system is represented at 
any instant by a state in a linear normed vector space, with a time-independent norm. Denoting by \i[>) a 
vector and its conjúgate partner by (ijjr¡\, the inner product is given by the volume integral, J é*r\édV = 
constant. Accordingly, the linearity of evolution and a normed space of functions implies 

d f 

— / é* v édV = 0, 
at J 

dé* f „ dé 



dt IY J r ' dt 

iip*H^r)tpdV - / ié*r]HédV = 0, 



J itp*r](r]- 1/ H t r]- / H)édV = 0, 

•H f = rflirf x . (1) 
We have made use of linearity (via, e.g., Schrodinger equation) 

Hi¡> = ihj£. (2) 

Let us consider symmetry transformations which preserve the 77-norm ((x\r¡y)) between the vectors x 

rH 

and y. By considering the Cayley form, D = e as a symmetry transformation acting on x, y where H is 
pseudo-Hermitian in accordance with r¡Hr¡ 1 = W, it is easy to show that D will be pseudo-unitary. Also 
it preserves the 77-norm and a consistently-defined matrix element provided D A D 1 = A', and is indeed 
a symmetry transformation ||9]. In general, metric for H and D are different and we will denote metric 
for H by 77 and for D by £. The metric will always be Hermitian but not unique. Additionally if metric is 
positive definite, then so is pseudo-norm. In contrast, any other metric leads to an indefinite pseudo-norm. 
The closure law for two pseudo-unitary matrices are trivial if they both are pseudo-unitary with respect 
to the same metric i.e. 77 = £. Also, D _1 is pseudo-unitary with respect to 77 if D is pseudo-unitary: 

r¡ (e )' r¡ = é iri i = e . With identity matrix as unit element of the symmetry transformation, 
and associativity guaranteed, the N x N pseudo-unitary matrices form a pseudo-unitary group of order N, 
PU(N). 

A pseudo-Hermitian operator possesses either real eigenvalues or the eigenvalues occur as complex- 
conjugate pairs. While this is trae in general, the essential point is regarding the corresponding eigen- 
functions. As discussed above, the right eigenfunctions have a corresponding left-partner which is related 
via 77 and conjugation. Precisely, for 'PT-symmetric system, complex conjugation corresponds to time- 
reversal and parity corresponds to 77. Connections of 2 x 2 matrices with 'PT-symmetry has been known 
1 9 . 52 . 56 1 . The eigenfunctions corresponding to the real (complex conjúgate) eigenvalues are (not) simul- 
taneously eigenfunctions of VT ■ The eigenfunctions corresponding to the complex conjúgate eigenvalues 
have a zero pseudo-norm [5|. Whereas the real eigenvalues correspond to bound states, the complex con- 
júgate pairs correspond to exceptional points. The latter are very interesting objects which have found way 
into our understanding of quantum phase transitions 11231 . zero-width resonances 11241 . quantum chaos Í25l . 
multichannel resonance ionization [261 etc. 

To understand the deviation from Hermiticity while still possessing real eigenvalues, let us take the 
example of 2 x 2 matrices where the eigenvalues can be written as 

1 „ . 1 



E± = -TtH ± - v/(TrH) 2 - ADetH. (3) 
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Now consider the example matrix, 



a —ice 
íc/e b 



(4) 



It is clearly not Hermitian as off-diagonal elements of the matrix are written as differently scaled elements 
of corresponding Hermitian matrix with e = 1. Now as we have not changed the trace and determinant of 
the matrix, from Eq. [3] the eigenvalues will be same to that of e = 1 (Hermitian) case and so they are real 
despite the matrix being non-Hermitian. 

Starting with the simplest case of a pseudo-Hermitian matrix IÍTÜI . 



H = {Hy} 



a —ib 
ic a 



a, b, c being real. It is easy to show that metric is 




i 

-i 



(5) 



(6) 



Interpreting this metric as parity (V), and usual complex conjugation, ICo as time-reversal operator T, the 
matrix is VT- symmetric. The diagonalizing matrix for H is given by D, Le., 



1 i/r 
ir 1 



(7) 



The eigenvalues of H are E± = a ± + -y] (r = \fc¡b (r G [0, oo] provided be > 0)). Here, as we 
have mentioned earlier that H(D) may not be pseudo-Hermitian (unitary) with respect to same metric, the 
metric for D is different from to that of H and is given by 



c = 



i 

1 o 



In the parameter space, the joint probability density function for the matrix H Í271 

P(H) = MeT 1 ^ 1 írHH 
is taken of the Wishart form, and this reduces to 



(8) 



(9) 



P(a,b,c) = ^ e -i^i 2a+h+c J. (10) 

V J 2(7TCT 2 )Í 

By inserting the relation of eigenvalues in terms of parameters and, using the Jacobian for this transforma- 
tion, we get the joint probability distribution function (jpdf) of eigenvalues: 

. \E+-E\ f(E + -E-) 2 \ (£++£^2l 
P(E + ,E_)= L± l Ko( K + - 2 M e 5==—. (11) 

Perhaps historically one of the most studied quantity in random matrix literature is the nearest neigh- 
bour level spacing distribution, P{S). It is well known that for the Wigner-Dyson ensembles the spacing 
distribution very well-approximated by, P{S) ~ 5^e _7S where ¡3 is 1, 2, and 4 corresponds to the or- 
thogonal, unitary, and symplectic ensembles [27, 28 29|. However, there are systems such as billiards in 
polygonal enclosures, three-dimensional Anderson model at the metal-insulator transition point, and many 
more those display intermedíate statistics ||3Ü1 13TI l32l . 
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The spacing distribution for the present case , P(S), is given in terms of the jpdf by 



/CSO /"OO 
/ P(E+,E-)6(S - \E+ -E-\)dE+dE- 
-oo J — oo 



151 K ( ^ 



(12) 



Remarkably, as S — > 0, P(S) ~ 51og(l/5) (it is non-algebraic level repulsión, in contrast with Wigner- 
Dyson ensembles). In TableJT] we summarize the results for various 2x2 pseudo-Hermitian matrices 
which form their own class. 

Table 1 A summary of 2 x 2 pseudo-hermitian random matrices 



H 



P(S) (S -> 0) 



a 


-ib 




' i' 


i 


"1 


i ' 

V 


ic 


a 




-i 


V2 


ir 


1 



(0 < r < oo) 



1 

1 



-SlogS 



a + c ib 




'1 " 


ib a — c 




-1 



eos 9 
—i sin í 



i sin# 
eos 6 



Veos 28 

(-7T/4 < 9 < 71-/4) 



Not known 



-SlogS 



-líC 

b 



a + ib c 
d a — ib 



e 

1' 

1 



coso lesmi 
-¿siné'/e eos 6 





\\ ol 




e 



sin í 

i 



re 

sin í? 
1 



Sf(i),e = e> 
Not known 5 with large slope 



a + b 


d- 


- ic 




"1 


" 




i eos 6 


é e sin 0" 




"1 


" 


—d + ¿c 


a 


- b 







-1 




e~ ie sin 9 


—i cosí? 







-1 



It is clear that unlike Hermitian cases, simply because of the presence of a larger number of parameters, 
a wide spectrum of spacing distributions are expected. The natural extensión of these results to the gen- 
eral PU(N) case remains open. However, we have exact results for a special kind of pseudo-Hermitian 
matrices, viz. eyelic matrices or circulants. 



3 Cyclic Matrices and random matrix theory 

Let us consider an N x N cyclic matrix with real elements, {a¿}: 



M 



ai a 2 
a2 a 3 



a N 
a/v-i 

"i 



(13) 
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It is important to note that this matrix is, in fact, pseudo-Hermitian (pseudo-orthogonal) with respect to rj 



that is, 



1 





1 





1 

1 





(14) 



M ' 



r)Mr) 1 . 



(15) 



Since r¡ 2 = identity, I, consistent with the earlier discussion, r) may be called "generalized parity". We can 
distribute the matrix elements using (7), P(M ~ cxp[— ylTr(M^M)]), and obtain an ensemble of random 
cyclic matrices (RCM) that is pseudo-orthogonally invariant. It is a well-known 041 that a cyclic matrix is 
diagonalized by a Fourier matrix, U, which is unitary: 



^ = ^exp^(i-l)(¿-l). 
The eigenvalues of M are given by [34- 1 



N 



Ei = a p exp 
p=i 



27TZ . 

— (p- 1)(Z - 1) 



(16) 



(17) 



(l = 1, 2, N), there is always one real eigenvalue as a sum of the elements, due to permutation symme- 
try. 

Let's start with the analysis for an ensemble of 3 x 3 cyclic matrices. To derive the joint probability 
distribution function (JPDF)of eigenvalues, as all the correlations are related to it we again need to calcúlate 
the Jacobian of transformation from parameters of matrices to the eigenvalues. We immediately see that 
tr M f M = 3 (o? + a 2 + aj). The JPDF of eigenvalues P({Ei}) can be written as 



P(E U E 2 ,EZ) 



-A(Ef+2\E 2 \ 2 



(18) 



As there are also two complex eigenvalues present, the definition of spacing is taken in the sense of 
Euclidean distance between the eigenvalues. With this definition of spacing, we may define S 23 ■= 
\E 2 - E 3 \ = V3(a 3 - a 2 ) as well as S* 12 := \E X - E 2 \ = ||(a 2 + 03) + '^-{a 2 - o 3 )|. Obviously, 
5*12 = as E 2 and E$ are complex conjúgate. Spacing distribution for the complex conjúgate pair, 
^0(^23) is given by 



r 3 

Pcc(S 23 ) = n[da i P({a i })S(S 23 -V3\ a 3 - a 2 \) 



(19) 



i=l 

Í2A 



A e¿ 

e 2 23 . 



The normalized spacing distribution by setting average spacing as 1, can be written in terms of the variable 

z: 



2 z 2 
Pcc(z) = -e - . 

TT 



(20) 
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Similarly, the spacing distribution, P rc {Si2) in normalized form can be obtained gain in in the variable 



z: 



3V37T 2 / 37T 2 2 \ /37T 2 o 

Prc(z) = c zexp c z io — c z 

y v ' 16 V 16 J V 32 



(21) 



wherec= 2 F x [f , |, 1, i]. 

It is clear that the Gaussianity of p cc (z) implies that there is neither a level repulsión ñor any attraction 
among the complex conjúgate pairs while real and complex eigenvalues display linear level repulsión. The 
numerical simulations are also in agreement as can be seen in Fig. [TJand Fig. [2] 





0.7 




0.6 




0.5 


-ti 


0.4 


o 




o 
CL 


0.3 




0.2 




0.1 










• numerics (3 x3) 

analytic 

x numerics (100 x 100) 







Fig. 1 Probability distribution of the absolute spacing between the complex conjúgate pair of eigenvalues of a Gaus- 
sian ensemble of 3 x 3 cyclic matrices. The numerical result matches well with the analytic result ( 120b . 



1 




12 3 4 

z 

Fig. 2 Probability distribution of the absolute spacing between a real and a complex eigenvalue of a Gaussian ensem- 
ble of 3 x 3 cyclic matrices (red stars). The numerical result obtained by considering 10000 realizations for 3x3 
matrices and 1000 realizations of 100 x 100 matrices (blue circles) agrees with the analytic result (solid black)(Eq. 
I2U . A difference from Wigner distribution for GOE (black dotted) is clearly shown in the figure. 
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In the same manner the jpdf for the general case of N x N matrices (with even N) is found to be 



A 

— | exp 

7T 



N 



A| E\ + E 2 n_ +í + 2^ EíEn+2-í 
i/1,4 +1 



(22) 



where E\ and En_ +1 real and the rest of the eigenvalues may be complex. For odd 7Y, there will be only 
one real eigenvalue, E\ and the summation in the second term will extend over all i except 1 . 

For the general case of 7Y, one more spacing will appear which was absent in 3 x 3 case and that is 
between two complex eigenvalues which are not complex conjúgate. Rest are found to be the same as in 
3x3 case. This spacing distribution which we denote be simply p(s) in normalized form is found to be 



p(s) 



y exp 



7TS 



(23) 



which is exactly the Rayleigh distribution [15| (Fig. |3). In contrast to Wigner's result which was an 
(excellent) approximation for N dimensional symmetric matrices, here it is an exact result for all N. 
Interestingly the same distribution occurs inthe case of a random Poisson point process in a plañe. Due to 
the fact that Wigner's distribution is approximate for N x N matrices, the above result is interpreted more 
gainfully in terms of the Rayleigh distribution of complex eigenvalues. 



100 x 100 
Wigner distribution 




Fig. 3 Probability distribution for spacing between two complex eigenvalues which are not complex conjúgate. The 
agreement with analytical result (Eq. I23) is good. 

As it has been already established that M is an example of pseudo-Hermitian (orthogonal) matrix, and 
as these results are found for N x N matrices, we believe that the results on cyclic matrices not only 
extend the random matrix theory in a significant way but also provide an illustrative example for more 
general results found for the Ginibre orthogonal ensemble with Gaussian distributed real elements solved 
recently by Í35l [36l . This also serves as an example where degree of freedom for the matrices are very 
constrained (only N), and how that affects the general resullts. 

4 Biased random walks on a regular lattice and cyclic matrices 

Let us first consider a random walk on a one-dimensional lattice of TV equally spaced sites with periodic 
boundary conditions. Assuming that the jump probability is w, it decides to jump to left orright neighbour 
with probabilities pw and qw(= (1 — p)w) respectively. Let us consider an ensemble of such lattices, and 
define the probability of occupation of the site, i by 



Pi = Ni/N e . 



(24) 
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where Ni denotes the number of lattices (realizations) with a partióle occupying the i site and N, 



N 



5> 



At time í, a state of an ensemble can be written as a vector 

i T 



Pi(í),...,Pjv(í) 



where T stands for transpose. The time evolution of the ensemble is given by 

p(t + 1) = Mp(t) 

where 



M 



(25) 



(26) 



(1 - w) 


pw 





qw 


qw 


(1 - w) 


pw 








qw 


(1-w) . 


pw 


pw 





qw 





(27) 



Matrix M is a transition matrix which can be easily recognized as an asymmetric cyclic matrix. Since this 
matrix is not Hermitian, its eigenvalues occur in complex conjúgate pairs, in addition to some of them being 
real. As we know, the Fourier matrix is the diagonalizing matrix for cyclic matrices. The n th component 
of the l th eigenvector corresponding to the eigenvalue, A; is 



(é¡)r 



1 



exp 



2ni 



(n-l)(l-l) 



(28) 



Now, any initial distribution can be expanded in terms of the eigenvectors of cyclic matrix and time evolu- 
tion simplifies to essentially taking powers of eigenvalues and then recombining which yields 



N 



(29) 



This equation is used to calcúlate the time evolution of Pi(t) for lattices of several different numbers of 
sites and for given valúes of w. 

We will use Boltzmann's relation to calcúlate the entropy of the system, S = ks ln fí with the thermo- 
dynamic probability, 



n 



N l 

1 > fJn .c ■ 



{Ni\)(N 2 \)...(Nm)' 



(30) 



Using Stirling's approximation for large N, and scaling by N ens , the ensemble averaged entropy, s 

S/N ens is 



s 



= - ^2 pí ín Pi 



(31) 



As equilibrium distribution corresponds to the largest eigenvalue, p¿ = 1/N s u e , and henee the limiting 
valué of s/ks is ln N s ue- The evolution of entropy is shown in Figj4] 
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100 



Time 

Fig. 4 For the case of random walkers on a periodic lattice with 22 sites, ai = 0.2, a-2 = 0.24, 022 = 0.56 (other 
üí's are zero) with the jump probability, w = 0.8, the entropy is seen here. 



5 Biased random walks on a disordered lattice 



We can generalize the previous model by allowing finite jump probability to all the sites, and let us in- 
troduce randomness in the elements of the transition matrix. The question is if transition matrix of the 
biased random walk is chosen from a Gaussian ensemble of cyclic matrices, then what can be said about 
the spectral properties and henee the evolution of entropy of this system? 

It is important to note that the matrix elements are probabilities and their sum is unity. Let us define the 
average by 



RCM = 



J¿drJ^d0p{t)p(r,O)r 
¡¿drf^d9p(r,6)r 



(32) 



where p{r,9) is the density of eigenvalues (A = re 10 ). It can be shown that r (= |A|) is distributed by 
a Wigner distribution (Eq. l23l l of random matrix theory, by the similar arguments as for Eq. |231 (also 
seeED). 

Using Eqf29]and Eqf28] pj (t) can be re-written as 



^¿>n(0) (AÍ + ¿Afexp 

n=l \ 1=2 

N N 



2-kí 



(l-l)U-n) 



(33) 



1=2 



where we have exp 



N 



= ilj n and also made use of = 1- Now, we can average 



over the the joint probability distribution function of eigenvalues, J(\2, Ajv)- Denoting the average as 
(Pj(t) - j¡{{^í}))rmt, we can write 
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N N 

(Pí(*)({Aí}))juít - ^E^WE^ 1 / d\ 2 ...d\ l ...d\ N J(\ 2 ,...,\ N )\ 

r~ — 1 í — n " 



n=l 
JV 



= 2 



i JV JV r i /•+(» 

r»=l ¡=2 L • y_oc 

¿V 



d\X t p(X) 



AT(AT- 



r* exp [if0] |r 2 exp (~ - ¿(é 1 ) 



kdrdO 



1 



1 



VE 

2 



JV ^_ e -7r/4 + Erf 
Using the identity for large a (V= 0, — 1, —2, ...) and fixed z, J40 

oo 

j{a,z) = z a e~ z Y^ 



2 ' 4 



z*T(a) 



r(a + fc + 1) 



and after a little algebra, we can rewrite N(pj (t) ({ A¿ }))hmt f° r large time í, 

2 7T 



JV(Pí(í)({A.-})) 



-e-^/ 4 + Erf 



2 



í + 3 (í + 3)(í + 5) \t 3 



(34) 



35) 




Fig. 5 With time p becomes more and more uniform, for three different N valúes the curves show how fast the 
uniformity sets in with time. 

The exclusión of first eigenvalue which always corresponds to = leads to the exclusión of the 
density of this angle (6(8)) from the uniform density of 6, valid for the rest of the eigenvalues. 

The time dependence of (Pj(t)({Xi})) rmt is shown in FigfS] Approach of all the probabilities to 1/N 
is evident, and so is the approach to a non-equilibrium steady state with máximum entropy given by log N. 
This is not surprising. Utilizing the connection cyclic matrix models, we believe that not only Gaussian but 
other kind of randomness can also be modelled. This a clean and simple way to do what could otherwise 
be done by setting up a master equation with a non-Hermitian Hamiltonian 1411 and solve for the steady 
state solution to study the approach to equilibrium. 
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(1.5 ■ 
Í0.4 
Z0.3- 

ta 




Bq.21 

Asymplolics i'or Eq. 21 



Fig. 6 After absorbing N in {Pj{t)({Xi})} rmt , the exact time dependence and first two term of asymptotic series 
expansión is plotted, in inset the difference between the two in % is plotted. It is clear that the convergence of the 
series is very fast. Asymptotically, the time dependence is given by < 



t+3 



6 Cyclic blocks, Ising Model and pseudo-Hermiticity 



In the two-dimensional Ising model defined on a regular square lattice, we associate a spin variable cr, with 
valúes +1 or -1 with each site. The nearest neighbour interaction terms is J^uíU^, and zero otherwise. 
With fixed J¿j = J, the partition function was found by Onsager [37, 38 ]. However, as is well-known (p. 7 
of ll27l ). if Jij is a random variable, with a symmetric distribution around zero mean , we have the random 
Ising model. the calculation of the partition function of which is an open problem. 

In the treatment of Onsager and Kaufman, a very crucial role is played by the transfer matrix with 
elements containing the coupling which was identified as having cyclic block structure. It would be nice to 
study the matrix models with such cyclic block structures. We have already seen that cyclic matrices with 
scaler entries are pseudo-Hermitian too. Let's see what can be said about cyclic blocks Í39l ? 

In this Section, we present a special case of random matrices with block entries. Firstly, let us consider 



B 



Ai A 2 
A N A 1 



A 2 A 3 
where each entry A¿ is 

Ai = 



A N 
A N -i 

Ai 



(36) 



-bi 



(37) 



o¿, bi, and c¿ are drawn independently from a Gaussian distribution with zero mean and unit variance. The 
matrix B is pseudo-orthogonal with respect to the "generalized parity", 



E = 



a 





a 



a ... 
where a is the Pauli matrix, 



(38) 



a 



1 

1 



(39) 
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Thatis, Bt = EBXT 1 . 

However, numerically it is seen that the spectral fluctuations are the same as that for random cyclic 
matrices with scaler entiles. For 50 x 50 matrix, comprised of fifty 2x2 blocks per row similar spacing 
distributions among complex conjúgate pairs, real-complex pair, complex-complex pair is found to have 
in good agreement with same in case of cyclic matrices with scaler entry i.e. Eq. |2Ü1 12T1 and |231 This is 
evident in Figs|7] and[8]respectively. 

10 ° . , , , 10° « , ^ , 



númenes (50 x 50 ) 
block of (2 x 2) 

- analytic 



numeric (50 x 50) 
block of (2 x 2) 
Eq. 11 from Ref. 
(Jain & Srivastava) 



Fig. 7 (Left panel) Log-log plot of the distribution of spacing among complex-conjugate pairs. The result is exactly 
in agreement with Eq. |2Ü]for scalar entries.(Right panel) Log-log plot of the distribution of spacing among real and 
complex eigenvalues. The result is exactly in agreement with Eq. |2T](Eq 11 from Ref. Jain & Srivastava) for scalar 
entries. 




Fig. 8 Log-log plot of the distribution of spacing among complex eigenvalues. The result is exactly in agreement with 
Eq. |23]for scalar entries (Eq. |2T]is Eq. 1 1 from Ref. Jain & Srivastava). 



The agreement is more intriguing than satisfactory as one would expect the block structure to show up 
in some form or the other in spectral distributions too. This led us to explore the cyclic blocks of (2 x 2) 
matrices having complex entries, numerically. In particular, keeping Ising model transition matrix in mind 
we focussed on cyclic matrices with the row (AB B . ..St) where 



.4 



ai ici2 
-ÍQ2 cti 



(40) 



B 



-I íbi 
ib 2 -r 



(41) 
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with real ai, 0,2, b\, 62. We have chosen the form of A and B from the stracture appearing in the Ising 
model in two dimensions. The spacing distributions for different cases for random matrices are shown in 
Figs E] [TÜ1 The complex-conjugate pairs are spaced as shown in Fig. |9]which is in reasonable agreement 
with Eq. [TU 




Fig. 9 Log-log plot of the distribution of spacing among complex-conjugate pairs. 

The spacing distribution between real and complex eigenvalue shows a departure in the tail (Fig. 
[TÜlhhowever we are yet to investígate the dependence of this deviation with the size of the matrix. Similar 
deviation of the numerical result from Wigner distribution is seen in Fig. [TÜ] 




Fig. 10 (left) Log-log plot of the distribution of spacing among real and complex eigenvalues. (right) Log-log plot of 
the distribution of spacing among complex eigenvalues (Eq. [2TJis Eq. 1 1 from Ref . Jain & Srivastava). 

This exploration throws a lot of interesting questions like obtaining the joint probability distribution 
function for the cyclic block models, and the ensuing spectral properties. 

7 Remarks 

Finally, we would like to point out a few significant áreas we have not touched upon. There is a lot of re- 
search directed towards the physics beyond the standard model, particularly relevant to atomic and nuclear 
physics. In line with the theme of the Special Issue is atomic electric dipole moment, the discovery of which 



Copyright line will be provided by the publisher 



14 



SCL Srivastava and SR Jain: pH-RMT 



will be an instance of breakdown of parity and time-reversal symmetries. In certain nuclei with shapes that 
would support such a symmetry-scenario, it is expected that certain collective modes would enhance the 
Schiff moment [42]. The quadrupole and octupole collective modes could co-exist. This área would gain 
if many-body theory is systematically extended beyond mean-field and random phase approximations, to 
include VT— symmetry in the semiclassical descriptions. The description of nuclear dynamics requires 
a self-consistent coupling of Hamiltonian for nucleons and a Hamiltonian in the deformation space. This 
leads to rather non-trivial connection of chaos in single-particle motion, emergence of hydrodynamic mode 
in deformation space giving rise to an effective inertia tensor l43l . 

Generally, there are much lesser known results on non-hermitian random matrix ensembles. The basic 
result is called the circular law, the analogue of the Wigner's semi-circular law for the density of eigen- 
values. This can be stated as spectral measure for n x n non-Hermitian iid matrices, whose elements are 
distributed in iid manner with fixed distribution (say Gaussian) with mean and variance 1, converges to 
uniform measure in circle. A figurative versión will be eigenvalues if scaled properly of such a matrix will 
fill the unit circle uniformly. For Gaussian unitary non-Hermitian ensemble the joint probability density 
and correlation function were obtained by Ginibre [44 1 while density of Ginibre orthogonal ensemble was 
given by Sommers et al. ¡45] and latter joint probability density function as well as the correlation function 
by Kanzieper and Akemann ll35l . A nice review on random matrices cióse to Hermitian and unitary is by 
Fyodorov and Sommers 11461 . In future, the assumption of joint independence of the matrix elements needs 
to be relaxed. Moreover, the non-Hermitian case presents the situation of spectral instability wherein a 
small perturbation in a large non-Hermitian matrix leads to large fluctuations. Similarly, there are a num- 
ber of interesting problems that remain open. For some new results on RMT for non-hermitian systems, 
we refer the reader to the work by Tao [l4Tl . 

The results on random 'PT-symmetric systems are expected to be relevant for a wide variety of physical 
situations occurring in anyon physics [48], i/-parametrized quantum chromodynamics, fractional quantum 
hall systems l49l . etc. The connection of RMT with the momentum distribution functions was established 
for the Maxwell-Boltzmann, Fermi-Dirac, and Bose-Einstein cases [50, 51 ]. As pointed out in [51 52 1, the 
momentum distribution function for the anyon gas requires the development of general pseudo-Hermitian 
RMT. Some current development related with T'T-symmetric deformations of integrable models is re- 
viewed by Fring lÍ53l while another interesting arricie on supersymmetric many-particle quantum systems 
with inverse-square interactions are now available ll5"4l . 

Non-hermitian circular billiard is shown to possess energy levéis that show some level repulsión, along- 
with the eigenfunctions which are directional 11551 . Further work on these lines will be relevant for the 
design of micro-disk lasers. 
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